Possible realization of the Exciton Bose Liquid phase in a hard-core boson model with 

ring-only exchange interactions 
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We investigate a hard-core boson model with ring-only exchanges on a square lattice, where a K\ 
term acts on lxl plaquettes and a K2 term acts on 1x2 and 2x1 plaquettes, with a goal of realizing 
a novel Exciton Bose Liquid (EBL) phase first proposed by Paramekanti et al [Phys. Rev. B 66, 
054526 (2002)]. We construct Jastrow-type variational wave functions for the EBL, study their 
formal properties, and then use them as seeds for a projective Quantum Monte Carlo study. Using 
Green's Function Monte Carlo, we obtain an unbiased phase diagram which at half-filling reveals 
CDW for small K2, valence bond solid for intermediate K 2, and possibly for large K2 the EBL 
phase. Away from half-filling, the EBL phase is present for intermediate K2 and remains stable for 
a range of densities below 1/2 before phase separation occurs at lower densities. 

PACS numbers: 71.10.Pm, 75.10.Jm, 75.40.Mg 



I. INTRODUCTION 

The Exciton Bose Liquid (EBL) theory by 
Paramekanti et al proposed a critical bosonic phase 
which shows remarkable resemblance to electrons in 
a metali For this novel quantum phase, the presence 
of loci of "Bose surface" in the Brillouin zone closely 
parallels the Fermi surface for fermions, and allows the 
Bose system to share many characteristics normally 
associated with fermions^ Paramekanti et al showed 
that due to the gapless lines of excitations, EBL is a 
critical (power law) compressible 2D quantum phase 
with uncondensed bosons and contains continuously 
varying exponents. Their striking proposal stimulated 
a number of works seeking to establish the stability of 
the EBL phase in bosonic models with ring exchange 
interactions^— However, these studies found that the 
EBL is not realized in the hard-core boson model on 
the square lattice with ring exchanges on elementary 
plaquettes. Instead, such ring interactions favor a (%, n) 
charge density wave (CDW) in the half-filled case, while 
away from half-filling they induce strong tendencies to 
phase separation. 

Motivated by the unusual EBL properties, our study 
focuses on a search for this unconventional quantum 
phase in simple bosonic models. A candidate model re- 
quires additional interactions for suppressing the charge 



( ~\ 

O 1 « 



-o 



• 




') p* c 




• 


— 


— « 


1 (b) 4 


I 


O 



G • (c) 



O 



FIG. 1: The ring exchange interactions for (a) lxl, (b) 2x1 
and (c) 1x2 plaquettes, which hop two bosons on opposite 
corners of a plaquette onto the two remaining vacant corners. 



order. One choice that comes to mind might be to intro- 
duce second nearest neighbor repulsion. In this paper, we 
instead adopt a different route where we consider addi- 
tional ring exchanges that frustrate the CDW tendencies 
of the elementary ring exchanges. We define plaquette 
exchange operators 
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where b r annihilates a boson on a site r, and x, y are the 
unit vectors on the square lattice. The Hamiltonian is 



h = —K\ J2 pi 1 - K 2 J2 ( p r 12 + p? 1 : 



(2) 



Figure Q] illustrates the action of these ring exchanges on 
hoppable plaquettes. The original ring model proposed in 
Ref.G] and studied numerically in Refs. 0-0 corresponds 
to K2 = 0. To sec how the present K1-K2 model may sta- 
bilize the EBL phase, we first note that the (7r,7r) CDW 
in the Xi-only model results from having a large num- 
ber of basis states connected to the perfect (7r,7r) CDW 
configuration. However, the K2 terms would be com- 
pletely inoperative in such a CDW. Furthermore, the P 12 
and P 21 ring exchanges by themselves would favor differ- 
ent charge orderings and not compatible with each other. 
Thus, the K2 ring terms compete with the K\ terms and 
with each other, making the liquid phase with no charge 
order more likely. 

The present K\-K<i model has the same lattice sym- 
metries and boson number conservation on each row and 
column as the original ring model of Ref. [H From the 
outset, we define our Hilbert space as the sector with 
equal number of bosons on each row and column. Note 
that this restriction does not preclude phase separation 
(PS); in fact, we shall see that PS does occur at low 
densities within our restricted Hilbert space. For non- 
negative K\ and K2 values, the Hamiltonian does not 
have a sign problem and allows an unbiased study of the 
system using Quantum Monte Carlo methods. Although 
the Stochastic Series Expansion is the method of choice 
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for simulating large lattices, it has not been applied to 
ring-only hard-core boson models due to implementation 
issues.— We instead use the Green's Function Monte Carlo 
(GFMC) approach with full bias control as described in 
Ref. 0, which allows us to obtain exact ground state prop- 
erties for moderately sized systems up to 12 x 12 in this 
work. Without loss of generality, we set K\ = 1 and vary 
Ki > in the study. 

Our main results for the phases of the model arc sum- 
marized in Figs. [TT] and [T4j First, at half-filling, our 
intuition that the K-i should suppress the charge order 
is indeed borne out, and the CDW disappears already 
for moderate K%. Somewhat surprisingly, this does not 
stabilize the EBL right away but instead drives the sys- 
tem into a columnar Valence Bond Solid (VBS), while 
the EBL is tentatively stabilized only for quite large K2 
terms. On the other hand, away from half-filling, mod- 
erate K2 already produce stable EBL phase. 

The paper is organized as follows. In Sec. UH we 
first construct good trial wave functions for the EBL 
phase and study their formal properties using Variational 
Monte Carlo (VMC), followed by an energetics study to 
determine the variational phase diagram. Using the op- 
timal trial states as starting point for the GFMC pro- 
jection in Sec. IIIH we compute density, plaquette, and 
bond structure factors in the ground states, followed by 
finite size scaling to determine the phases. Our phase 
diagram reveals the (tt,tt) CDW order at small K2, VBS 
order for intermediate K2, and possibly for large K% the 
novel EBL phase. We perform detailed comparison of 
the numerical results with the EBL theory. In Sec. IIV1 
we extend the search for the EBL to densities less than 
half at intermediate -K2, and find that the EBL liquid is 
stable for 1/3 < p < 1/2 while phase separation occurs at 
lower densities. In Sec. [Vj we conclude with a discussion 
of possible further studies. For the benefit of readers, 
Appendix summarizes the results of the EBL theory 
from Refs. Ill relevant for our numerical work. In 
Appendix [B] we offer a parton-gauge perspective of the 
EBL theory. 



II. VARIATIONAL STUDY AT p = 1/2 
A. Formal properties of the EBL wave function 

In this section, we study the formal properties of the 
EBL wave function. To provide motivation for the wave 
function, we first consider a quantum rotor version of the 
model with elementary lxl ring exchanges, 
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where the phase 4> r and the boson number n r are canoni- 
cally conjugate. In the EBL theory, the cosine in the ring 




FIG. 2: Density structure factor Ssw(qx,q y ), Eq. for the 
ground state of the spin- wave Hamiltonian, Eq. (j4j). The char- 
acteristic "cross" formed by the singular lines q x — or q y = 
is a distinctive feature of the EBL phase. 



term is expanded to quadratic order (this approximation 
is valid in the stable "spin-wave phase" with no topolog- 
ical defects) . The resulting set of coupled harmonic os- 
cillators can be diagonalizcd in momentum space, which 
leads to the following, 
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We will loosely refer to ifsw as the "spin-wave" Hamil- 
tonian. Writing its ground state in the n variables and 
then restricting to n r £ {0,1}, we obtain a valid hard- 
core boson wave function in the convenient Jastrow-type 
form that can be implemented easily in VMC, 
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In the spin wave theory, W = \/U /K, while here it serves 
as a variational parameter. For the K\-K% model, W 
becomes a q-dependent function with two parameters. In 
this section, we will focus on the single-parameter EBL 
wave function to illustrate properties of such variational 
states and what can happen with them. 

To characterize the phases realized in the EBL wave 
function, we measure the density structure factor 
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The density structure factor of the ground state for the 
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spin-wave Hamiltonian in Eq. (|4]) is given by 
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At any fixed q y , S$w(q x ,q y ) vanishes for small q x as 
C{q y )\q x \, with further C(q v ) ~ \q y \ as (fy — > 0. This 
gives the characteristic "cross" shown in Fig. [5] which is 
a signature of the compressibility of the EBL, and is re- 
lated to the gaplessness of ujq along the loci q x = or 
q y = 0. The latter is a consequence of the conservation 
of boson number along each row and each column of the 
lattice. To identify possible realization of the EBL phase, 
we monitor the long wavelength behavior of the density 
structure factor in addition to the absence of Bragg peaks 
in all structure factor measurements made in this paper. 

For sufficiently strong interaction U (in particular, for 
hard-core model) at half-filling, the spin-wave approxi- 
mation no longer holds. Proliferation of topological de- 
fects results in the (tt,tt) CDW instability found in earlier 
Quantum Monte Carlo studiesi^r— Remarkably, we find 
that the single-parameter EBL wave function Eq. ^ is 
able to realize both the EBL and the CDW phase. As we 
increase W in the half-filled system, the wave function 
undergoes a phase transition at a critical value W c ps 4.4 
where the (n, it) charge order develops. This is analyzed 
in the top panel of Fig. [3] using finite size scaling of 
S(tt,tt)/L 2 , which vanishes as 1/L 2 in the absence of the 
order for W < W c and approaches a finite value in the 
presence of the order for W > W c . The non-monotonic 
L dependence of this CDW order parameter for fixed 
W > W c is somewhat unusual but appears to be a prop- 
erty of such wave functions, perhaps indicative of some 
long crossovers in the system. 

Next we examine the long wavelength behavior of the 
density structure factor near the characteristic "cross" . 
Figure [3] also shows the ratio 



v{Qx,Q v ) 



S(q x ,q y ) 
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evaluated at the smallest q x = q m i n = 2tt/L as a func- 
tion of q y on lattices with length L between 10 and 80. 
For W = 4.2 in the middle panel, the ratio shows some 
deviation from csw = 1/(2 W) but it clearly renormal- 
izes towards finite values as expected in the EBL theory. 
When plotted in the full Brillouin zone, the VMC density 
structure factor looks essentially like Fig. [5] 

On the other hand, for W = 4.6 in the bottom panel, a 
strong downward renormalization of the ratio is observed 
for all q y , in particular near q y = 0. Such behavior is simi- 
lar to a Mott insulator, where the density structure factor 
is nonsingular and hence has cuts (q x — > 0,q y ) with van- 
ishing slopes. Thus, the contrasting long-wavelength be- 
haviors of S(q Xl q y ) independently confirm a phase tran- 
sition in the wave function near W c ~ 4.4. 

To summarize, the above wave function with one varia- 
tional parameter can realize either the EBL liquid phase 
or the (tt, it) charge order and thus can alert us about 
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FIG. 3: (color online). VMC study of the EBL wave function 
in Eq. © with one parameter W defined by Eq. ([TJ- Top: 
Finite size scaling of the density structure factor S(7r,7r). The 
wave function undergoes a transition at a critical value W c ~ 
4.4, which separates the EBL phase at low W from the (tt,tt) 
CDW at larger W. Middle and bottom: a(q x =2ir/L,q v ), 
Eq. ()10|) . which gives normalized slopes of the density struc- 
ture factor near the cross, plotted against q y for lattice length 
L = 10 to 80. Middle: Results for W = 4.2 show long wave- 
length EBL characteristics where these normalized slopes ap- 
proach finite values [apparently close to crsw = 1/(2W) indi- 
cated with a dotted line]. Bottom: Results for W = 4.6 show 
a downward renormalization of the slopes, which is similar 
to a gapped state. The critical values estimated from both 
the order parameter study and the detailed study of the cross 
agree, i.e., W c « 4.4. 
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FIG. 4: Finite size scaling of the VMC density structure factor 
S(n, 7r) obtained in the energetics study using optimized W 
for each K2 in the Hamiltonian and each system size L. 



CDW tendencies in the system. A note of caution is ap- 
propriate here. Our GFMC study in later Sec. IIIII shows 
that the formal wave function study of the present sec- 
tion does not always capture a full physics of the problem. 
Namely, as we will discuss in Sec. IIII C 11 the wave func- 
tion itself may be in the liquid phase, while the full EBL 
theory with the same effective K/U is already unstable. 
Nevertheless, our formal wave function study clearly has 
its own merits. For instance, it alerts us to the possibil- 
ity of complex crossovers with the system size and that 
the order may be weak and not apparent on short scales, 
but may still appear on longer scales. It also teaches 
us to look at the long-wavelength behavior for signs of 
instabilities. 
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FIG. 5: Variational phase diagram for the K1-K2 model on 
the half-filled square lattice, based on the single-parameter 
trial wave function. Refer to Fig. [11] for a revised phase dia- 
gram based on unbiased GFMC results. 



III. UNBIASED GFMC STUDY AT p = 1/2 

In this section, we perform a Green's function Monte 
Carlo study which, being unbiased, provides an impor- 
tant check on the VMC results. In the GFMC approach, 
an initial trial wave function is projected onto the ground 
state via repeated application of a projector which even- 
tually filters out the excited states: 



|V^+i) = [1-(H- E )St] 



(11) 



Here Eq is a parameter chosen close to the ground state 
energy, and St is a "time step" chosen to ensure domi- 
nance of the ground state and positiveness of the projec- 
tor, which then allows Monte Carlo calculations without 
a sign problem. Operator expectation values are eval- 
uated using stochastically sampled ground states which 
generally requires the so-called GFMC "forward walk- 
ing" technique, and we implement this using the bias- 
controlled scheme described in Ref.0- We refer the reader 
to the literature for more details on the GFMCiii~— 

To identify the nature of the true ground states, we 
measure the density structure factor S(q x ,q y ) defined in 
Eq. © as well as the following plaquette structure factor 



P(qx,q y ) 



1 



,iq-(r- 



i(p T "y (p r vn , (12) 



B. K1-K2 energetics study with one-parameter 
EBL wave function 



Let us turn to the energetics study using the above 
one-parameter wave function. For each value of K 2 and 
system size L we find the optimal variational parameter 
W (later Fig. |5] compares the trial energies with the ex- 
act diagonalization values on the 6x6 lattice) . In Fig. 21 
finite size scaling of the corresponding structure factor 
S(n, 7r) shows the presence of the CDW order for small 
K 2 and the absence for K 2 > 0.3. Thus, the VMC study 
suggests that the EBL phase could be stabilized even 
with quite weak K 2 ring interactions. Allowing an addi- 
tional variational parameter corresponding to K 2 in the 
spin-wave Hamiltonian does not modify this conclusion. 
Figure [5] summarizes the VMC results obtained for the 
K\-K 2 ring model at half-filling. In the following section, 
our GFMC simulation reveals another phase in the inter- 
mediate K 2 region. We will present a revised phase dia- 
gram based on unbiased GFMC results in Sec. IIII B[ and 
will discuss the failure of the VMC later in Sec. IIII C II 



where (-P,. 11 ) 2 equals 1 if the lxl plaquette is "hoppable" 
and otherwise. While quantitatively different from the 
off-diagonal P} 1 plaquette structure factor used in Rcf. 3, 
the operator (P,. 11 ) 2 defined here is easier to implement in 
the GFMC and it gives qualitatively the same access to 
bond-solid-type phases. To better discriminate between 
plaquette and bond orders, we also measure the following 
bond structure factor 
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where B? = btb r 
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is 1 if the bond is "hoppable" and otherwise, and is 
again easy to implement in the GFMC. 



A. Test of our GFMC setup 

In Fig. [BJ we test our GFMC setup against exact di- 
agonalization (ED) calculations for the 6x6 lattice. All 
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FIG. 6: Top: Comparison of the 6x6 lattice VMC and GFMC 
energies per site against the ED ground state and first excited 
state values for < K 2 < 0.9. For K 2 > 0.2, the VMC 
energies are already closer to the ground state than to the 
first excited state; however, despite such a good agreement 
in the energies, the VMC fails to identify a new phase for 
K% > 0.3 as described in the text. Bottom: Comparison of 
the corresponding VMC and GFMC structure factors S(n, n) 
against the ED values. The GFMC energies and structure 
factors show complete agreement with the ED values. 



results are in the sector with 3 bosons in each row and 
in each column. On both panels, the GFMC and ED re- 
sults essentially coincide for the entire K2 range shown. 
Table [J summarizes the respective data for reference. To 
check the accuracy of the trial wave functions, the VMC 
energies are also plotted in the top panel of Fig. [6] along 
with the first excited state ED energies. For K2 > 0.2, 
the VMC energies are already closer to the ground state 
than to the first excited state. We caution the reader 
that despite this very good accuracy on the 6x6 system, 
the VMC fails to identify another order that develops for 
larger K2 and is found by the GFMC for larger sizes. 
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TABLE I: Comparison of the ground state energy and S(n, tt) 
obtained using the VMC, GFMC, and ED calculations for the 
6x6 lattice. The energies are given in units of K\ per lattice 
site. The GFMC results are essentially exact, and we treat 
them as such for larger sizes. 



B. GFMC study of the K1-K2 model at half-filling 

We now proceed to apply this numerical tool to char- 
acterize the ground states of the K\-K<i ring model. The 
top panels in Figs. [7] and [8] show the density structure 
factor £(71", 7r) and the plaquette structure factor P(7r,0) 
plotted against K% for lattice sizes ranging from L = 6 
to 12. Between K2 = and 0.4, S(tt,tt) increases 
strongly with L while the size dependence weakens with 
K2 ■ This coincides with a strengthening size dependence 
of P(7t,0). Beyond K2 = 0.4, the charge order is absent 
while the plaquette order now dominates in the range 
up to K2 ~ 4. For still larger K2, P(tt,0) becomes 
very weakly dependent on lattice size. We do not ob- 
serve any other strong feature in S(q x ,q y ) and P(q x ,qy) 
over the full Brillouin zone. Thus we identify the (/T,7r) 
CDW for < K 2 < 0.4, a (tt, 0) bond-solid-type phase 
for the intermediate K2 region, and tentatively an EBL 
phase for K2 > 4. The middle panels in Figs. [7] and 
show finite-size scalings of the respective order parame- 
ters which support these conclusions. 

The bottom panels in Figs. [7] and [8] show the corre- 
sponding Binder ratios defined as 



Binder ratio 



(\M\ 4 



(\M\*)* 

j(7T,7r)-r 



(14) 



with Meow = Er e%{n,nyr nr or MvBS 
£ r e i(7r < 0) - r (P r n ) 2 (so (|Af| 2 ) are simply propor- 
tional to the already discussed structure factors). The 
Binder ratios provide additional information about 
the fluctuations of the order parameters (via the mea- 
surement (|Af| 4 )) and are expected to approach 1 in 
the presence of the order and 3 in the absence of the 
order. Such change in the behavior is clearly seen 
when the CDW order disappears and the plaquette 
order appears near K2 ~ 0.3 — 0.4. Note, however, 
that the familiar Binder crossing technique apparently 
does not work for the CDW order parameter for our 
sizes, even though we arc confident that the CDW order 
disappears (also supported by the values w 3 of the 
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FIG. 7: Top: GFMC density structure factor S(ir, ir) versus 
K2 for periodic lattices with length L = 6,8, 10 and 12. Note 
that to show the data compactly, taken for K2 = to 1 in 
steps of 0.1 and K2 — 1 to 10 in steps of 1, we used linear 
scale for the first range but log scale for the second range. 
Middle: Finite size scaling of S(tt,tt). Bottom: Binder ratio, 
Eq. <|14[) . for the CDW order parameter. Note the apparently 
large finite size effect in the ratio, particularly for sizes L — 6 
and 8 at large K2, and also absence of clear Binder crossings. 
Nevertheless, the Binder data is generally consistent with a 
lack of CDW order for K2 ^ 0.3, as can be seen from the ratio 
approaching the expected "disordered" value 3 (shown with 
a dotted line). 



FIG. 8: Top: GFMC plaquette structure factor P(n, 0) versus 
K2 for periodic lattices with length L = 6,8, 10 and 12. Linear 
scale is used for K2 — to 1 in steps of 0.1, and log scale for 
K2 — 1 to 10 in steps of 1. Middle: Finite size scaling of 
P(-7r,0). Bottom: Binder ratio, Eq. (|14p . for the bond-solid 
order parameter. The Binder data is consistent with no VBS 
order for small K2 and VBS order for intermediate K2. The 
observation of Binder ratios exceeding the disordered value 3 
suggests that there is no order for large K2, although there 
is no clear crossing going to the large K2 phase (note that 
similar somewhat unusual behavior is also seen in the CDW 
Binder data when the CDW order disappeared.) 
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Binder ratio itself) . Note also strong and non-systematic 
size dependence, particularly for the smallest L = 6. 
In the plaquettc Binder ratio, we see lack of order for 
small K%, appearance of order for intermediate K2, and 
apparently "disordered" Binder values for K2 > 4 (which 
is consistent with the absence of the plaquette order), 
but no clear crossings for this transition. Although the 
Binder data does not clearly give us the critical value 
of K2 for the transition to the disordered phase, its 
limiting value strongly suggests that there is no (tt, 0) or 
(0, 7r) plaquette order at large K 2 - An additional lesson 
from this study is that we should be aware of particular 
strong finite size effects in this system. 

Our identification of the VBS order for the intermedi- 
ate K2 region is further helped by measurement of the 
bond structure factor B x (q x , q y ) defined in Eq. ([T3]) . For 
this measurement (not shown), similarly to the plaque- 
tte structure factor, we observe (tt, 0) and (0, tt) Bragg 
peaks but no peak at (tt, tt). This is more consistent with 
a "columnar" VBS order rather than a plaquette order, 
and is also similar to the phase found by Sandvik et al£ 
in the J-K model for 8 < K/J < 14. Our finding of the 
same VBS state may in fact be related, since J added 
to the pure K\ model may induce effective K2 ring ex- 
changes frustrating the CDW while still remaining in the 
Mott insulator. 

To get a more complete picture, we examine the long 
wavelength behavior of the density structure factor us- 
ing the "cross analysis" of Sec. Ill Al Figure [5] shows the 
"normalized slopes" a(q x =2-K j L, q y ), Eq. (TT0|) . for L = 6 
to 12. The left panel shows the results for K2 = 0, which 
we already know is in the CDW phase from the presence 
of the (tt,tt) Bragg peak. We clearly observe a Mott-like 
incompressible behavior where the slopes vanish. This 
is similar to the earlier formal wave function study with 
the CDW. The Mott-like dependence of S(q x ,q y ) at long 
wavelengths continues to be present after the charge or- 
der disappears for K2 ^ 0.4. This is illustrated in the 
middle panel for K2 = 1, at which the bond-solid order- 
ing is already established in Fig. [5] Our "cross analysis" 
therefore provides an independent detection of the insta- 
bility to a different solid. 

In the right panel of Fig.|H]for K2 = 7, our small- lattice 
data appears to suggest that S(q x ,q y ) does have the V- 
shaped singularity along the lines q x — or q y ~ 0. This 
would mean that the bond-solid ordering exists only at 
intermediate K2 and hence, possibly realizing the EBL 
phase at large K2- However, this panel also reveals a 
weak downwards renormalization of the slopes upon in- 
creasing L and we therefore do not rule out the possibility 
of the EBL behavior disappearing at much larger lattice 
sizes. 

A rigorous confrontation of the large K2 region requires 
a study on much larger lattices, but this is beyond the 
capability of our present numerical setup. We instead ex- 
amine the structure factors S(q x ,q y ) and P(q x ,q y ) over 
the entire Brillouin zone and look for signatures of the 
EBL phase as well as potential instabilities. The top 



panel in Fig. [10] shows the density structure factor at 
K 2 = 7 and clear absence of any CDW ordering. Here, 
we highlight the presence of the long wavelength EBL 
signature near the lines q x = and q y = (this charac- 
teristic cross has already been analyzed in Fig. [9J. The 
middle panel in Fig. |TU] shows the plaquette structure 
factor for the same system, which again does not show 
bond or plaquette ordering. Despite the potential insta- 
bility hinted by the P(7r,0) and P(0,7r) cusps, the size 
independence of the plaquette structure factor along the 
cut q x = tt shown in the bottom panel of Fig. [TO] gives 
us some confidence that the EBL phase may indeed be 
realized in the large K2 regime at half-filling. 

Figure [TT] summarizes the unbiased GFMC phase dia- 
gram obtained for the K1-K2 model at half-filling. 



C. More detailed comparison with the EBL theory 
and interpretations 

We now discuss how the presented results at half- 
filling may fit into the EBL theory framework reviewed 
in Appendix [A] The EBL is characterized by an "EBL 
phase stiffness." For example, in the spin wave the- 
ory Eq. (|4]), the EBL stiffness is simply K/U, which 
we will parametrize by k = \J K/U. More generally, 
the EBL stiffness is a function on the cross lines q x = 
or q y = 0, with convenient paramctrization by n(q y ) = 
-\//C(0, q y )/U(0, q y ) in the notation from Appcndix[A] As 
described there, n(q y ) determines exponents in various 
power-law correlations along lattice directions, which in 
turn determine stability of the EBL. 

We do not have a direct access to the EBL phase stiff- 
ness in our setup, but we can crudely monitor its behavior 
using the characteristic cross in the density structure fac- 
tor. Specifically, from the EBL theory result Eq. (|A6|) . 
we have for the "normalized slopes" Eq. ([T0|) : 

^(0,%) = ^x|C p (0,g,)| 2 . (15) 

Here C p (q) is a non- universal function of order 1. In the 
spin-wave theory, C p (q) = 1. In general, C p (q — > 0) = 1, 
and we conjecture also that C p (q x 0,q y ) = 1 for any 
q y , but we do not know for sure. If we can indeed ignore 
the |C P (0, q y )\ 2 factor in Eq. (|T~5]) . we can then view our 
"cross analysis" presented earlier as a crude measure of 
the effective EBL stiffness on the corresponding length 
scales. When the measured a(q m \ n ~2ir / L , q y ) renormal- 
izes down to small values upon increasing the system size, 
the EBL is unstable; when a(2ir / L,q y ) stays finite, the 
EBL is stable. 

The theoretical stability of the EBL requires the stiff- 
ness to be sufficiently large, and the condition is partic- 
ularly stringent at half-filling because of allowed Umk- 
lapp interactions. The corresponding scaling dimensions 
are given in terms of n(q y ) in Appendix [A] For a q y - 
independent k, the leading Umklapp has scaling dimen- 




FIG. 9: "Cross analysis" of the density structure factors, plotting "normalized slopes" a(q x = 2n/L,q y ), Eq. (|10[1 . versus q y 
for L = 6 to 12. Left: K2 = in the CDW phase; the vanishing of the slopes is consistent with nonsingular behavior expected 
in an incompressible system. Middle: K2 = 1 in the VBS phase; the vanishing of the slopes continues to be present and can be 
used as an indication of the EBL instability even if we did not know the resulting order. Right: K2 = 7 (note different vertical 
scale); the data shows non-zero normalized slopes and only weak renormalizations, therefore suggesting a stable EBL phase. 



sion given by Eq. (|A"23|) and is irrelevant if 

K > 3/8 , (16) 

which we can take as a rough guide at half-filling. Our 
EBL cross analysis of the density structure factor gives 
a w k/2, so to establish stability we would like to see 
a > 3/16 = 0.1875. For K 2 = 7 this is satisfied on 
average for sizes L = 6 to 12, but the larger sizes come 
close to the threshold. The state may be somewhat more 
stable in that the measured cr(q y ) is largest near q y = n/2 
and the particular "average" of n(q y ) that one needs has 
the main weight in the middle of the region [0,7r], cf. 
Eq. (|A21[) . [We also want to repeat that we are not sure 
whether we can ignore the factor \C P \ 2 in Eq. ([To]) and 
unambiguously extract n from a.] From the absence of 
any orders, we conjecture that this point is stable; of 
course, if the instability is very weak, we may be not 
detecting the order on our length scales. 

For K2 = and K2 = 1, the effective EBL stiffness 
parameters as estimated by the cross analysis in Fig. 
are below critical or border-line critical already for the 
size L = 6 and then quickly renormalize down upon in- 
creasing L, consistent with our finding of the instability 
of the EBL towards boson solid phases. 

As reviewed in more technical terms in Appendix 
when the EBL is unstable at half-filling, the natural out- 
comes are a (n, it) CDW or a plaquette solid with period 
2 in both lattice directions. The outcome depends on 
the sign of some effective couplings. If an effective near- 
est neighbor repulsion dominates, the solid locks into the 
CDW, while if the second-neighbor repulsion dominates, 
the solid locks into the plaquette state. We then propose 
that, as we increase the K2, while the EBL remains un- 
stable (e.g., as detected by the cross analysis), at some 
point the sign of the locking switches from the CDW to 
the bond-solid. (At present, we do not know how to real- 
ize the columnar VBS out of the EBL theory, but usually 
columnar and plaquette orders are related^— and per- 
haps we are missing some physics ingredients in the the- 



ory that would enable the columnar VBS.) As we further 
increase K2, we conjecture that the Umklapp eventually 
becomes irrelevant and the stable EBL is realized. In this 
scenario, we do not anticipate any other instability, so if 
the presented large K 2 region is eventually unstable, the 
simplest possibility is that it will have a very small VBS 
order. 



1. Interpretation of the failure of the VMC at p = 1/2 

In light of the above stability considerations, we now 
briefly discuss the failure of the formal wave function 
study in Sec. [TT| to detect the EBL instability in the in- 
termediate K 2 regime. We presented mainly the one- 
parameter wave function that can capture only the EBL 
or CDW. However, we also considered a two-parameter 
wave function where in the spin wave theory like Eq. (0| 
we include both K\ and K2 ring terms (in fact, this was 
used throughout to obtain improved initial states for the 
GFMC projection). We found that such a wave function, 
depending on the parameters, can realize also the VBS 
state on the same footing as the CDWJ^ Nevertheless, 
in the energetics study in the intermediate K2 regime, 
the optimized two-parameter wave function produces a 
liquid. 

The one-parameter example from Sec. [TT] is sufficient 
for our discussion. The wave function is constructed from 
the EBL spin wave theory and it seems reasonable to take 
the EBL parameter as n = 1/W. Consider now W = 4.2 
shown in the middle panel in Fig. [31 where the normal- 
ized slopes in the density structure factor approached the 
expected value crsw = 1/(2 W) and where we concluded 
that the VMC wave function is in the liquid phase. How- 
ever, such stiffness K = 0.24 strongly violates the stability 
condition Eq. (|16p . Therefore, we appear to have a sit- 
uation where the formal wave function is a liquid with a 
stiffness that is too small for the full EBL theory to be 
stable. 
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FIG. 10: Top: GFMC density structure factor S(q x ,q v ) for 
a 12x12 lattice at K 2 = 7. Middle: GFMC plaquette struc- 
ture factor P(qx,q y ) for the same system (q = point not 
calculated). Bottom: Cut of P(q x ,q y ) at q x = n for L — 6 to 
12. 
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FIG. 11: GFMC phase diagram for the K\-K 2 model on the 
half-filled square lattice. We found the (tt, 7t) CDW for K 2 ^ 
0.4, the (n, 0) VBS for intermediate values of K 2 , and possibly 
the EBL phase for K 2 > 4. 



This is reminiscent of what happens when one formally 
considers a Jastrow-type wave function for ID hard-core 
bosons, if! id = IIi<j I sin[7r(xj — Xj)/L]\ v . The wave func- 
tion describes a Luttinger liquid of bosons with the Lut- 
tinger parameter g = \jv. On a half-filled chain, the 
Luttinger liquid becomes unstable to a staggered CDW 
when g < 1/2 corresponding to v > 2. However, as dis- 



cussed in Ref. Il9l the above wave function remains liquid 
until v exceeds 4, and only then the CDW order develops. 
Thus, in this ID example, the condition for the formal 
stability of the wave function is different from that in the 
full theory. 

Assuming similar phenomenon for the formal EBL 
wave function, we can then speculate what happened in 
our variational study. The wave function parameter is 
found by optimizing the energetics and it roughly cap- 
tures the bare EBL stiffness on the scale of few lattice 
spacings. For small K 2 , this is already in the regime 
where both the wave function and the full EBL theory 
are unstable. However, for intermediate K2 the opti- 
mized parameters happen to be in the range where the 
wave function is stable while the full EBL theory is not, 
hence the failure of our VMC. Of course, once we sus- 
pect boson-solid phases, the use of the few-parameter 
wave functions motivated from the liquid side becomes 
inadequate. In the variational approach, more param- 
eters also allowing the Jastrow pseudopotentials to be- 
come more long-ranged would be needed ; 20 ' 21 while in 
the present study the correct physics is brought by the 
GFMC projection. 



IV. STUDY OF THE K x -K 2 MODEL FOR p < 1/2 

When we step away from half-filling, the Umklapp 
terms discussed above are no longer allowed. While the 
EBL may still be unstable due to non-Umklapps, they 
are typically less relevant (see Appendix EJ. However, 
here one also competes against phase separation at low 
densities. Previous studies^^ of boson models with 
lxl ring exchanges found that ring interactions induce 
strong tendency to phase separation, since they are oper- 
ative only when bosons are nearby. The more extended 
K2 ring interactions can somewhat offset this tendency 
and produce a stable uniform EBL regime over a wider 
range of densities below half-filling. 

We would like to point out that our restricted Hilbert 
space with equal boson number in each row and each col- 
umn does not preclude phase separation. For example, 
basis states with preferential clumping along a diagonal 
or in blocks along the diagonal are present in our Hilbert 
space. In fact, we observe regimes of phase separation 
in the VMC and GFMC simulations which will be dis- 
cussed below. We detect the phase separation in the 
Monte Carlo simulations either by monitoring snapshots 
of the real-space boson configurations, or by looking at 
the structure factors in momentum space, where it is re- 
vealed by the presence of strong peaks at the smallest 
wavevectors. To further check the results, we start our 
simulations from both uniform (random) boson configu- 
rations and from half-filled diagonal stripes. We verified 
that, independent of the initial configurations, our sim- 
ulations converge to uniform states for densities close to 
1/2 but phase-separate for low densities. 

It may be true that in our working sector the phase scp- 
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Pbos 

FIG. 12: Trial energies of the optimized EBL wave function 
measured on a 24x24 lattice. Using Maxwell construction 
(illustrated for K2 — 1), the critical density dividing the uni- 
form liquid regime and the phase-separated regime can be 
determined for each K2. We found a decrease in the critical 
density from 0.45 to 0.37 as K2 increases from to 1. 

aration is somewhat suppressed in finite samples, since 
other shapes of clumped regions are not allowed. How- 
ever, this effect should decrease with increasing system 
size, and it is likely that our sizes already crudely capture 
such local energetics tendencies as to whether the system 
wants to stay uniform or phase-separate. 

Below, we present results of a VMC energetics study on 
a 24x24 lattice and results of a GFMC study on lattices 
with L < 12. We find that a stable EBL phase is present 
in a window p <£ (0.4,0.5) for K2 as small as 0.5. 

A. VMC results for 24x24 lattice 

For our VMC energetics study of the uniform liquid 
phase on the 24x24 lattice, we use the single-parameter 
EBL wave function from Eq. ©. We consider total bo- 
son numbers in multiples of 24, as appropriate for the 
Hilbcrt space with equal number of bosons in each row 
and column. Figure [12] shows the trial energy per site 
for the optimized EBL wave function plotted against bo- 
son density. Away from half-filling, uniform VMC liquids 
are obtained for all data points shown in the figure. For 
p < 0.13, we already detect phase separation in the VMC 
simulations (data points for these cases were excluded 
from the plot). Among the uniform liquids shown on the 
plot, not every point produces a stable liquid. We can use 
Maxwell construction (illustrated for K2 = 1.0) to deter- 
mine the critical density p c dividing the stable uniform 
liquid regime (p c < p < 1/2) and the phase-separated 
regime (p < p c ), the latter comprising a uniform liquid 
region and an empty region on the lattice. We found 
a decrease in the critical density p c from approximately 
0.45 to 0.37 as K2 increases from to 1.0. Thus, the 
additional K2 ring terms indeed help to widen the stable 



regime of the uniform EBL phase in the K1-K2 model. 

We should of course be cautious taking the VMC re- 
sults too literally, given the described experience with 
the failures of the VMC at half-filling. However, we are 
probably in a better position here in that the "bare EBL 
stiffnesses" are such that both the wave function and the 
full theory are stable. Specifically, for K2 > 0.2 and all 
densities p < 1/2, the optimal W is smaller than 3.5 and 
is further decreasing with increasing A2. In the absence 
of Umklapps, the most important non-Umklapp has scal- 
ing dimension 8y/K/U = 8/W > 2 (sec Appendix I A 3|) ; 
hence, all residual interactions are irrelevant and the full 
EBL theory is stable. While the VMC results are sug- 
gestive, the ultimate determination of the phase diagram 
requires unbiased approaches. 



B. GFMC results for L < 12 lattices 

Next, we discuss the results of our GFMC study on the 
12 x 12 lattice for two selected points in parameter space, 
K2 = 0.5 and 1.0. We consider densities varying from 2 
up to 5 bosons per row and per column (corresponding 
to densities p = 1/6 up to 5/12). For 2 and 3 bosons per 
row at K2 = 0.5 and 2 bosons per row at K2 = 1.0, we 
already see signatures of phase separation both in real 
space and momentum space. For higher boson number 
per row and per column, we do not see any sign of the 
phase separation. This allows us to conclude that the 
12x12 system is in a stable uniform phase (i.e. without 
phase separation) for 1/3 < p < 1/2. 

To determine whether the uniform phase realizes the 
EBL, we examine the GFMC density and plaquette struc- 
ture factors for any sign of instability to CDW or bond- 
solid ordering, but we do not observe any strong peak 
in S(q x ,q y ) or P(q x ,q y ). We apply the "cross analysis" 
of earlier sections to study the long wavelength behavior 
of the density structure factor. The top panel of Fig. [TBI 
shows the analysis done for K2 = 1-0 and p = 1/3, which 
is carried out on L x L lattices with L = 6, 9 and 12. 
Note that the crude stability condition at generic den- 
sities is a « k/2 = (l/2)y/ K/U > 1/8 and is safely 
satisfied. Note also that the bare EBL stiffness on the 
scale L = 6 is similar here and in the half-filled system 
with the same K 2 shown in the middle panel in Fig. |H1 
however, unlike the half-filled case, it does not renormal- 
ize upon increasing L consistent with the picture where 
some relevant Umklapp becomes inoperative for p < 1/2. 
(One needs to worry about higher-order Umklapps if the 
density happens to be commensurate, which we do not 
worry here, having more in mind incommensurate densi- 
ties in a window 1/3 < p < 1/2. Density p = 1/3 may 
be slightly outside the stability window if we take the 
VMC energy per site estimates in Fig. [T2] seriously and 
is perhaps stabilized here against phase separation by the 
finite system size, but is a good example allowing us to 
see the absence of flow of the EBL stiffness with several 
our sizes.) Similar result is obtained for A'2 = 0.5 (not 
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FIG. 14: GFMC phase diagram for the K\-Ki model on the 
12 x 12 lattice with varying boson density from 2 up to 6 
bosons per row and per column. 



which includes densities below 1/2 that we studied. We 
have also added points at K2 = 0, where our similar 
study using the "cross" technique indicates that the EBL 
is unstable also close to half-filling (but we have not es- 
tablished the resulting phases). Thus, it has been crucial 
to add moderate K2 exchanges to realize the EBL away 
from half-filling. 



CONCLUSIONS 



FIG. 13: The GFMC density structure factor S(q x ,q v ) on 
the 12x12 lattice at density p = 1/3 and K2 = 1.0. Top: 
"Cross analysis" showing normalized slopes a{q x = 2tv/L, q y ), 
Eq. I |10| l. versus q y . Bottom: S(q x ,q y ) in the full Brillouin 
zone. Besides the "cross" signature, notice also (weak) "2/cf 
ridges" [here lines (±27r/3, q y ) and (q x , ±2tt/3)] which are typ- 
ically present for all the uniform liquids observed at density 
less than half. 



shown), and both indicate that strong EBL signature is 
present. In the bottom panel of Fig. Q21 the sharper 
"cross" shows that the EBL is more stable for density 
p < 1/2 compared to the half- filled system. 

We also note the presence of small 2kp ridges (dis- 
cussed in Appendix which we typically observe in the 
GFMC density structure factor of uniform liquid ground 
states studied at densities away from half. This feature 
is predicted in the EBL theory and may be taken as ad- 
ditional evidence for identifying the uniform phase with 
EBL. Thus, our GFMC study for densities close to 1/2 
shows that the intermediate K2 regime is a stable EBL 
phase. Together with the VMC results, we think the ev- 
idences are sufficiently strong to conclude that the EBL 
phase is already realized in the K1-K2 model for densi- 
ties 1/3 < p < 1/2 and intermediate K2 values, while the 
phase separation dominates at lower density. 

Figure [TI] shows the extended GFMC phase diagram 



In summary, we studied the K1-K2 hard-core boson 
model with ring-only exchanges and found a transition at 
half-filling from a staggered CDW order at small K2 to 
a columnar VBS order at intermediate K2 , and tantaliz- 
ingly realizing the EBL phase at large K 2 - For densities 
away from half- filling (but not too far), the EBL phase 
is more robust and our evidence strongly suggests that 
the EBL is already realized in this model for intermedi- 
ate K2- (Having some K2 is helpful, since the pure K\ 
model does not appear to have the EBL also away from 
half- filling.) For still lower densities, instability to phase 
separation is observed. 

Although our sizes arc significantly smaller than typi- 
cally studied for more conventional phases with advanced 
techniques such as Stochastic Series Expansion, we al- 
ready reach 12 x 12 systems which are much larger than 
sizes used in ED studies. Like the ED, the GFMC 
method that we use provides exact information about 
the ground state and allows us to reasonably establish 
the phase diagram of the proposed model already with 
our sizes. Our GFMC evidence of the EBL phase is quite 
suggestive at half filling, although we critically point out 
possible pitfalls. Specifically, while direct measurements 
of the charge- or bond-solid order parameters do not 
reveal any orders, the detailed comparison of the long- 
wavelength signature in the density structure factor with 
the available EBL theory suggests that our tentative EBL 
points at half-filling are very close to being unstable. (It 
may be helpful to modify the model still a bit like adding 
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2x2 ring exchanges to more reliably stabilize the EBL in 
the half- filled system.) On the other hand, similar com- 
parison for the claimed EBL points away from half-filling 
indicates that they are safely away from instabilities. 

We hope that it will be possible to study significantly 
larger sizes and confront our tentative EBL realizations 
as well as systematically confront the EBL theory. There 
are a number of further properties and questions that one 
would like to explore. Thus, we have measured only den- 
sity and energy correlations that are diagonal in the num- 
ber basis as they are simplest to implement in the GFMC. 
The gaplessness of the EBL and its "Excitonic" charac- 
ter can be probed directly by measuring boson "box" 



At 



correlations like Gy(x,y) = (&( 0l o) 6 (o,tf)&(s lI /) & (*.o)) 
fixed y, this can be viewed as an exciton propagator for 
excitons of size y, and is predicted to show power-law be- 
havior ~ with calculable y-dependent exponent.— 
A crucial characterization of the EBL (which, in partic- 
ular, determines all power law exponents) is the "EBL 
stiffness" function. While we have had some access to it 
via the cross analysis of the density structure factor, it 
would also be interesting to measure the EBL stiffness 
directly. More broadly, the EBL is an example of a very 
gaplcss quantum liquid, and it can be challenging but 
fruitful learning grounds for how to handle such phases 
in Quantum Monte Carlo simulations. 

One immediate question for the EBL theory is that 
in its present form it does not seem to anticipate the 
columnar VBS phase that we found numerically in the 
K1-K2 model. Another question for both numerical and 
theoretical studies is to understand phases away from 
half- filling near K2 = 0, in the regime where we did not 
find the phase separation but also concluded that the 
EBL is not stable. 

We also note that our present realizations of the EBL 
are very likely immediately unstable towards a supcrfluid 
if we allow unfrustrated boson hopping. This is based on 
our estimates of the EBL phase stiffnesses and under- 
standing of the EBL stability conditions against boson 
hopping.— Thus, an EBL phase envisioned in Ref. [H in 
such a broader sense [with no special conservation laws 
other than global U(l)] will not be realized with our K\- 
K2 model. However, we hope that our work will stimu- 
late further studies to achieve such a phase. 

Our tentative numerical realization of the EBL (even 
in the restricted sense of ring-only systems) is of broader 
interest in the search for so-called Bose-metal phases 
or more generally non-Fermi liquid metals.— ~— Despite 
many studies (including more recently of holographic 
metals in the High Energy Theory communit y 31 ! 32 ), to 
date there is no example where such phases can be 
demonstrated controllably. The EBL can be viewed as a 
special kind of a Bose-metal and its theory is on more firm 
ground (an interesting perspective on the EBL is that it 
can be viewed as a solvable example of a parton-gauge 
theory where partons have flat Fermi surfaces in mean- 
field, cf. Appendix |B1) . It is hoped that the present work 
may trigger more refined studies of the EBL and con- 



frontation with the theory, and more efforts to access the 
challenging but very interesting and topical Bose-metal 
phases. 
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Appendix A: Precis of the EBL theory 

For the benefit of readers, this Appendix summarizes 
results from Rcfs. Hi relevant for our numerical work, 
complementing them where needed. 



1. Gaussian fixed point theory 

The EBL phase is described by a fixed point 
Lagrangiani 



Co = H [cp,#}- 
/C(k) 



^ -c^MA^R) , (Al) 

7T 



\{A X y(p)l 



M(k) 
2tt 2 
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We use the same notation for the fields as in Ref. [H ip 
corresponds to the spin wave part of the original phase 
variables, while d describes coarse-grained boson density 
fluctuations <5n(r) = n(r) — n. In the microscopic lattice 
derivation^ 

5n Q (r) = - |tf(R) - t?(R - x) - tf(R - y) (A2) 



■tf(R-x-y) 



= -A xy $ 

TT 



(A3) 



Lower-case r = (x, y) denote the original lattice sites, 
while upper-case R = (X, Y) denote the dual lat- 
tice sites, cf. Fig. fTS) The above Wo contains general 
"EBL parameters'— /C(k) and U(k), which are func- 
tions of k = (k x ,ky) respecting square lattice symme- 
tries. Long-wavelength properties such as various power 
law exponents (discussed below) depend on the func- 
tion JC(0, k y )/U(0, k y ). We will consider stability of this 
fixed point to allowed perturbations later. A parton- 
gauge derivation utilizing familiar ID Bosonization rules 
is sketched in Appendix [Bl 

Integrating out the field i9 gives a spin wave theory 
Lagrangian in terms of the if field that corresponds to 
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the spin wave Hamiltonian Eq. ((4]). However, the com- 
pactness of the original phase variables (or, equivalently, 
the discreteness of the boson number) is best captured in 
terms of the dual field Therefore we consider a dual 
Lagrangian obtained upon integrating out the tp field, 



C 



dual 



" (k) |(A^) k | 2 



1 



27T 2 IV XU /K| 27T 2 /C(k) 

(A4) 

Calculations mentioned below are performed in this 
Gaussian theory, and the main task is to express physical 
operators in terms of the continuum fields. 

The Gaussian theory gives equal-time correlations of 
the coarse-grained boson density as^ 



|<5n (k)| 




sin(fc x /2) sm(k y /2)\ (A5) 
(A6) 



The "cross" formed by singular lines (0, fc y ) and (k x ,0) 
is the most prominent feature in the density structure 
factor in the EBL phase. Note that the result contains a 
non- universal factor |C p (k)| 2 , which is expected to arise 
after coarse-graining 5n (r) — > J2 m C p (m)Sno(r + m). 
The coefficients C p (m) are non- universal but satisfy 
E m Cp(m) = 1, implying C p {k -> 0) = 1. We con- 
jecture that specification of the EBL fixed point theory 
(including identification of observables and residual in- 
teractions below) requires an even stronger condition: 
C p {k x — > 0,k y ) = 1 for any k y . 

A more accurate formula for the boson number op- 
erator contains also ll 2kp v contributions (where 2kp = 
27m), 

n(r) = n + 5n a + A cos[2V y i? + 2nn(x - 1 /2)] (A7) 
+ Acos[2V K z? + 27rn(y- 1/2)] , (A8) 



where n is the boson density per site and 
V M tf = i?(R) - i?(R - A) • 



(A9) 



[More precisely, the oscillatory contribution say in 
the first line in Eq. (| A7|) should be written as 
sin[2Vj,i?(X, y) + 2-xnX] - sin[2V y tf (A -l,y) + 2nn(X - 
1)], where the original lattice coordinate x is between 
dual lattice coordinates X — 1 and X as in Fig. [T3J How- 
ever, any fuzzincss in the location of along the x di- 
rection tuns out to be unimportant and the more precise 
expression simplifies to that in Eq. (|A7|) .] The partic- 
ular oscillatory contributions are not explicitly listed in 
Ref.[H but require a more sophisticated treatment there; 
they are discussed in Refs. 00 and are also recapped in 
Appendix [Bl from the parton perspective. 

To characterize the 2kp oscillatory part of the boson 
density, we need correlations 



/ e i2V M tf(R) e -i2V„tf(R/) 



>Y,Y> 



\ X _ X ,V2A 2 



(A10) 



These are ultra-short-range in the y direction. The power 
law in the x direction (as well as in the imaginary time 
direction not written explicitly) is determined by the scal- 
ing dimension 



= 2 



/K(0,M 
U%ky) 



sm{k v /2)=A\j-. (All) 



For an illustration, the very last equality specializes to 
the case of k-independent K. and IA ; here and below this 
case is marked with "=". Note that in the actual den- 
sity correlations, Syy' is replaced by an exponentially 
decaying function of \Y — Y'\, since by general symmetry 
reasoning we also expect contributions in Eq. (|A7I) where 
V y i9 is evaluated at r±jy, albeit with exponentially de- 
creasing amplitudes. In k-space, the density structure 
factor has line singularities, 



\k x ± 2nn 



2A-2 



(A12) 



In the parton perspective on the EBL, Appendix [Bl these 
can be viewed as the 2kp surfaces associated with partons 
moving in the x direction. There are also similar 2kp 



lines at k y = ±2irn. 



Finally, the same singular lines arc present in bond and 
plaquette correlators. For example, for an energy-type 
operator associated with a bond [r, r + x], we havei^ 



r,r+x 



A' cos[2V y t? + 2nnx] 



(A13) 



which is similar to the oscillatory part in the density op- 
erator except for a phase shift. 



2. Nonlinear interactions 

The discreteness of the original boson charges can be 
faithfully represented by allowing non-linear terms in the 
■d variables. The first such interactions that one would 
write down aroi 



C = - J2 °<? cos[g(2tf + 2nnXY)} 



(A14) 



Such insertions have ultra-short-range correlations when 
calculated in the Gaussian action Eq. (|A4|) . and small 
such residual interactions are therefore irrelevant in the 
EBL theory. However, as discussed below, there are al- 
lowed additional interactions (which can be viewed as de- 
rived from the C n ) that do have power law correlations, 
and then stability requires that the power laws are suffi- 
ciently fast. Note that non-integer boson density causes 
oscillatory phases, which can effectively "disallow" inter- 
actions from the low-energy theory. Note also that even 
though the terms in Eq. (|A14[) do not enter the stability 
condition for the EBL fixed point, once the EBL becomes 
unstable they are important to correctly describe result- 
ing phases (see analysis in Ref. [l], Sec. V B 2, and here 
in Sec. lA4a[) . 
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3. Stability of the EBL at generic densities 

We first discuss stability of the EBL phase at incom- 
mensurate density n. For any n, there is an allowed (non- 
oscillatory) interaction associated with chains parallel to 
x and separated by jy: 

A = -fi 0¥0) cos [2V„i?(R) - 2V a tf(R + jy)] . (A15) 

There are similar interactions associated with chains par- 
allel to y. Such interactions were not discussed in Ref. 0], 
while j = 1 case was considered in Ref. HI (see also Ap- 
pendix |EJ . The above operator has short-range correla- 
tions in the y direction, but power law correlations in the 
x direction as well as in the imaginary time; the corre- 
sponding scaling dimension is 



Aft 



Mi - 




rJS?44 sin 2 (jfc s /2) sin(fc y /2) 



(A16) 



The last line specializes to the case of k-indcpcndent fC 
and Li ; some representative values are 



A[h^ =±1) ] 

A [ f) b = ±2)] 



Interestingly, in this case the smallest scaling dimension 
is achieved for j — > oo, probably because of particular 
inter-chain correlations present in the system. 

Without any assumptions on the parameters /C(k) and 
Li(k), the scaling dimension in the j — > oo limit can be 
related to the 2k p correlation dimension; 




A[f)fr' =0 °)] = 2A 



2h r . 



(A20) 



References [l|HH argued that stability to such residual 
interactions is determined by comparing the scaling di- 
mensions to a reduced space-time dimensionality equal 
to 1 + 1 = 2 here. Adopting this criterion, we conclude 
that stability requires A 2 fc F > 1. Therefore, the 2kF 
correlations must decay faster than 1/x 2 along chains, 
or, equivalently, the corresponding singularity Eq. (|A12[) 
across the 2kp lines in the momentum space must be 
weaker than slope discontinuity. 

It is not clear at present what happens if the system is 
unstable to the above interactions in the incommensurate 
case, i.e., in the absence of Umklapp terms. Another 
remark is that if one also allows boson hopping between 
chains and considers k-independent K, and Li, then the 
scaling dimension for the boson hopping works oulji to be 
such that there is no stable EBL phase with unfrustratcd 



hopping even at incommensurate density. It would be 
interesting to explore different EBL stiffness functions 
^//C(0, k y )/Li(0, k y ) that may give stability against both 
Eq. (|A15[) . and boson hopping. In the present paper, 
we focus solely on the restricted class of models with no 
boson hopping. 



4. Stability of the EBL at half-filling 

Let us consider stability at density h = 1/2 where large 
part of our numerical work is performed. In this case, the 
following interactions are also allowed, 

C u = -u ij) cos [2Vj,n?(R) + 2V y tf(R + jy)} . (A21) 

Here j = corresponds to an Umklapp in each chain, 
while general j corresponds to an Umklapp involving 
chains separated by jy (cf. Appendix [Bj. There are simi- 
lar interactions associated with chains parallel to y. The 
corresponding scaling dimensions, in the same sense as 
before, are 



Afu (j) l = 



i 




(A22) 



The last line again specializes to the case of k- 
independent K, and Li. Here the smallest scaling dimen- 
sion is obtained for j = ±1 corresponding to interactions 
between nearest neighbor chains, 




A[u^=°>] = 16 W — . 

V Li 

The last two lines give the scaling dimensions of 
the second-neighbor-chain Umklapp and the intra-chain 
Umklapp. 

Reference [l| considered the following Umklapps at half- 
filling, 



-ro 2g , m cos [2g0(R) - 2qd(R + my)] , (A26) 



with the condition that qm is an even integer. We have 
tt>2,2 = u*- 1 ) and tn^i = u^°\ but in general the sets 
u' J ) and K)2q,m arc different. The latter Umklapps can 
be viewed as involving m neighboring chains (cf. Ap- 
pendix [Bj . The scaling dimensions arei 



A[ro 



2g,mJ 



2f 



4^ 



1 1 

1+ 3 + 5 



IJC(0,k v ) sm 2 (mk y /2) 
U{Q,k y ) sin{k y /2) 

' + 2m -1 J V U 
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where the last line specializes to the case of k- 
independent K, and IA. We also caution that more 
complex interactions not considered here may get re- 
duced scaling dimensions for some functional forms of 

K(0,ky)/U(0,ky). 

As far as boson ring models are concerned, the EBL 
fixed point is stable for sufficient dominance of the effec- 
tive "ringing" K, over the "repulsion" U. 

a. Instabilities at half-filling 

Let us consider the case when the most relevant term 
is the nearest-chain Umklapp tD2,2 = uS 1 ' [this is true, 
e.g., for k-independent K/U, see Eq. (|A23[) ]. We also 
assume u*- 1 ' > 0, which is expected in the presence of 
effective nearest-neighbor repulsion, see discussion after 

Eq. ram 

It is convenient to define 

9 1 {X,y) = V y d , (A27) 
9 2 {x,Y) = V x ti . (A28) 

In principle, all discussion below can be carried out using 
i? variables only and is essentially equivalent to the anal- 
ysis m Ref. [0, Sec. V B 2; one addition is discussion of 
interactions and regimes that can produce a pure CDW 
state. Nevertheless, the &\ and 2 variables introduced 
here make a connection with the parton-gauge perspec- 
tive on the EBL, Appendix [Bj and, up to a point, allow 
to use quasi-one-dimensional language and intuition. For 
example, the interactions Eqs. (|A15[) and (|A21[) are sim- 
ply expressed in terms of the 9\ or 62 and are naturally 
associated with parallel chains in one or the other direc- 
tion. However, note that these variables are constrained 

by v x 0i =v y e 2 . 

When the U^ 1 ' Umklapp is relevant and flows to strong 
coupling, it provides pinning 

20i (X, y) = -20 1 (X,y-l) = (-l)«20 1 (X,O) (A29) 
= (-1)^20! (0,0) , (A30) 
26 2 {x,Y) = (-1)* 20 2 (0,0) , (A31) 

where all equalities arc modulo 2tt, which is natural 
periodicity for the 2t9, 29\, and 202 variables. In the 
above, we have also minimized the quadratic energy 
Eq. (|A4[) by looking for time-independent fields satisfying 
V x 0i = V y 02 = 0. For the $ variable we obtain 

2&(X,Y) = 20(0,0) + ^— -20i (0,0) (A32) 

+ ^y— ^ 20 2 (O,O) . (A33) 

We are thus left with three undetermined parameters 
20(0,0), 20i (0,0), and 20 2 (O,O). 

We can loosely view the nearest-chain Umklapp uS 1 ' as 
locking relative density fluctuations in neighboring chains 



but not locking the density to the lattice. The latter can 
be provided by several terms which we list roughly in 
the order of their relevance in the EBL theory and which 
potentially have large bare values: second-chain Umk- 
lapp —vS 2 " 1 cos[20i(X, y) + 28 1 (X, y + 2)}; near-chain non- 
Umklapp -f)W CO s[20i(A, y) - 20i(X, y + 1)]; and intra- 
chain Umklapp — u^ ' cos[40i(A, y)]. We have only shown 
terms associated with chains running in the x direction; 
there are also similar terms associated with chains in the 
y direction. Putting everything together, we obtain an 
effective pinning potential on the 0i(O,O) and 02(0,0), 

-Aefr{cos[40i(O,O)] + cos[40 2 (O,O)]j . (A34) 

Schematically, ignoring any intermediate scale RG flows, 
A off =u( 2 )+b (1) +u(°> . (A35) 

Let us first consider the case A h > 0. This gives pin- 
ning 

20j (0,0) = tttoi , (A36) 

202 (0,0) = 7TTO 2 , (A37) 

where integers mi, m,2 can take independent values and 
1. By examining observables Eq. (|A13[) . we can already 
conclude that the system has columnar dimer patterns 



on the x- and y-bonds: 

&8r,r+x = yl'(-l) x cos[20i(X,y)] (A38) 

= A'(-l) x (-l) mi , (A39) 

SB r:r+ y = A'{-l)y{-l) m - . (A40) 



The simultaneous presence of these bond energy patterns 
corresponds to a plaqucttc solid. To complete the analy- 
sis, we also need to consider pinning of the 2$(0, 0). This 
is provided by terms in Eq. (|A14|) . For simplicity, con- 
sider just the Di term. Using Eqs. (|A33|) and (|A37|) . this 
energy is proportional to — t)i cos[2i?(0, 0)] [1 + (— l)™ 1 + 
(-1)" 12 - (-l)" ll+ " i2 ]. We readily establish that for each 
pair {mi,m 2 }, the 2$(0,0) is uniquely determined and 
gives a state where three out of four plaquettcs have more 
negative energies, i.e., a plaquette solid; four distinct 
pairs {mi, m,2} correspond to four ways to put this solid 
on the lattice. (Original Ref. H also considered adding D2 
to the energy mix and argued that for sufficiently neg- 
ative X>2 the system will have coexisting plaquette and 
CDW orders.) 

Let us now consider the case A c ff < 0. This gives pin- 
ning 

20i (0,0) = l+Trmj , (A41) 

20 2 (0,0) = |+7rm 2 . (A42) 

By examining observables Eq. (|A7[) . we can already con- 
clude that the system has (tt, it) CDW order: 

fei(r) = A(-irsm[20i(X,y)] (A43) 
= A(-l) x+y (-l) mi , (A44) 
Sn 2 (r) = A(-l) x+y (-l) m2 . (A45) 
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To be more precise, we also need to consider the pin- 
ning of th e 2^( 0,0) by terms in Eq. (|A14|> . Using 
Eqs. (|A33[) and (|A42j) . the Di energy is proportional to 
-Oi{cos[2^(0,0)][l + (-l) mi+m2 ]+shi[2i9(0,0)][(-l) mi + 
(-1)™ 2 ]}. This is minimized when m\ = m 2 and unique 
2-0(0,0) in each case, and produces uniform plaquette 
energies. The two independent minima mi = m 2 = 
or 1 correspond to two ways of putting the CDW on 
the lattice. [A note on approaches: Our Eq. (|A7[) as- 
sumes some coarse-graining, i.e., some high-energy fields 
are already integrated out. We could instead proceed 
using bare lattice variables $ where Sn(r) = A xy {}/TT 
and minimize the dual action including the nonlinear 
terms Ui, U2, f)i, Uo, 0i, etc. We would then find that a 
careful treatment of the saddle point approximated by 
Eq. (|A42[) indeed has a staggered charge density wave, 
5n(r) ~ (-l) x +y.) 

At this stage, we can speculate that the ring model at 
K 2 = has small bare value of u 2 and large bare values 
fli < 0,Uo < due to repulsive effective nearest- neighbor 
interactions, see discussion after Eq. (|B13I) . In this case, 
A c ff < and the CDW phase is realized. As we increase 
K 2 , it likely feeds directly into bare 112 > because of the 
effective second-neighbor avoidance desired by such ring 
terms. This can eventually make A c ff > and change 
the order to plaquette solid. However, our numerics sug- 
gests that the K\-K 2 model at intermediate K 2 realizes 
the columnar VBS instead. We have not succeeded to 
understand this within the above EBL instability treat- 
ments. It would be good to clarify this since the two 
bond-solids are usually related^— and perhaps we are 
missing some physics ingredients. [For example, one can 
contemplate more interactions, say leading to terms like 
-A'cos[40i(O,O)]cos[40 2 (O,O)] in addition to Eq. (|X34)| . 
which can indeed produce columnar VBS but only coex- 
isting with the CDW.] 



Appendix B: Parton-gauge perspective on the EBL 

Here we offer a parton-gauge perspective on the EBL 
phase. An effective "dimensional reduction" noted by 
previous authors^i^ 3 - 3 - - — can be also related to the 
one-dimensional (ID) character of partons. The par- 
tons are still strongly interacting, but here we can treat 
all gauge fluctuation effects accurately and in fact ar- 
rive at the EBL description summarized in Appendix |A"1 
We will also see how the familiar ID Bosonization 
technique a 36 i 37 allow to quickly obtain physical observ- 
ables and important residual interactions in the EBL 
theory. [We emphasize, however, that the EBL phase 
is qualitatively different from sliding or crossed-sliding 
Luttinger liquids 3 ^— — for example, it has specific heat 
C~Tln(l/T), cf. Rcf.CO.] 

We write 



(X-l.Y) (xjY) (X,Y) 

* * 



(X-l,y) (x,y) (X,y) 
□ O □ — 



(X-l.Y-l) (x, 



Y-l) (X.Y-1) 



n 



FIG. 15: The parton phase fields ipi(x,y),ifi2(x,y) reside on 
the sites of the original square lattice (black circles); the ID 
dual field 6i(X, y) resides on the horizontal links and 62 (x, Y) 
on the vertical links (white boxes). The dual EBL theory field 
&(X, Y) resides on the plaquettes of the original lattice, or 
equivalently on the sites of the dual lattice (stars). 



and recover the physical Hilbcrt space by the constraint 



n(r) = m(r) = n 2 (r) 



(B2) 



6t(r) = & t 1 (r)4(r) 



(Bl) 



We will arrive at the EBL theory by starting from a 
"mean field" where b\ partons hop only in the x direction 
while 62 hop only in the y direction^ One can justify 
such a starting point, e.g., by noting that this mean field 
gives large negative boson ring energies. As is familiar in 
slave particle treatments, a theory of fluctuations beyond 
the mean field contains a gauge field a x (r),a v (r) resid- 
ing on the links of the lattice; an alternative route con- 
necting the parton-gauge system and bosonic ring model 
can be found in Sec. Ill A of Ref. [U Here the partons 
bi and b 2 carry opposite gauge charges with respect to 
a = (a-xiO-y)- A complete theory also needs to treat the 
constraints Eq. (|B2[) . Typically this would be done by 
introducing an auxiliary field which would then be in- 
terpreted as a temporal component of the gauge field. 
However, here we are able to work with the constraints 
without the need for such a new field. 

We use phase and dual variables (pi(x,y) and 0i(A", y) 
to describe a harmonic flui d 36 ' 37 of b\ partons moving on 
a horizontally oriented chain located at a vertical coordi- 
nate y. Similarly, we use variables <p 2 (x,y) and 02 (a;, Y) 
to describe a harmonic fluid of b 2 partons moving on 
a vertically oriented chain at a horizontal coordinate x. 
Note that ipi(x,y) and ip 2 (x,y) reside on the sites of the 
original lattice, while 0\{X,y) resides on the horizontal 
links and 2 (x, Y) on the vertical links as illustrated in 
Fig. 1151 Including coupling of the partons to the gauge 
field, the Gaussian part of the imaginary time Lagrangian 



17 



reads 



£[v?iA, V?2, 2 ,a x , a y ] 



(B3) 



J2 [^( V ^2+a y ) 2 + 




d T ip 2 V y 9 2 



f ( v * a y - V y a K ) 2 



Here for simplicity we showed the nearest-neighbor "par- 
ton hopping" coupling J and on-site "repulsion" u, but 
all analysis below can be carried out more generally. Cru- 
cially, we assume a stable phase where (a x ,a y ) can be 
treated as a non-compact gauge field with a large "stiff- 
ness" parameter k (this assumes strong energetics selec- 
tion of the particular mean field state by the ring ex- 
changes of the microscopic boson model). The parton 
densities are 



nx(x,y) 
n 2 (x,y) 



V X X 



'y02 



so the constraint Eq. (|B2[) reads 



(B4) 
(B5) 

(B6) 



At this stage, we can integrate out the fields (pi and 
ip 2 and obtain Lagrangian density 

[[e u 6 2 ,a x ,a y ] = ^ [(V^i) 2 + (V y 2 ) 2 ] 

[(3 T 0i) 2 + (d T e 2 ) 2 ] 



2tt 2 J 



+ ~(a x d T 9i - a y d T 8 2 ) + -(V^ - V y a x y 

7T I 



Next, we solve the constraint Eq. (|B6|) via 

0i (X, y) = W y d = §{X, Y) - §{X, Y-l), (B7) 
9 2 (x, Y) = V x d = 0(X, Y) - 0(X - l,Y) , (B8) 

cf. Fig. |T5] Integrating out the field a, we finally obtain 
Lagrangian density 



u —rl „a\2 , 1 /a T7„cn2 



r(Vi.<?)" + 



2tt 2 J 



2tt 2 k 



(S T tf) 2 .(B9) 



This is essentially the EBL theory written in the dual 
variables Eq. Ki]) . with U(k) = 2u and 1//C(k) = 
4[sin 2 (fe r /2) + sin 2 (fc. i/ /2)]/ J + 1/k. [Note that one may 
be tempted to drop the (<9 T V$) 2 /J term as it contains 
more derivatives than the (9 t $) 2 /k term. However, the 
long-distance EBL properties such as power-law expo- 
nents are determined by the full function /C(0, k y ) which 
does depend on J if we want to be accurate in the simple 
model Eq. (|B3|) that we took.] If we include from the 



start general interactions among the partons and general 
Maxwell terms for the gauge field, we obtain the gen- 
eral Gaussian EBL theory with k-dependent /C(k), U(k) 
described in Appendix [Al 

We can now establish connections between microscopic 
observables and the EBL theory in the # variables. Thus, 
in the ID Bosonization treatment, the particle density 
and bond energy are given by 



5ni (r) 



A cos [20! +27rn(a;- 1/2)] (BIO) 



A' cos [26»i + 27rnx] 



(Bll) 



Here 9\ already means the long-wavelength component 
and the precise location where it is evaluated along the 
chain is unimportant. Writing 9\ via Eq. (|B7[) we obtain 
Eqs. (|A7[) and (| A13|) quoted in the EBL theory precis. 

We can also express inter-chain density-density inter- 
actions 

Vjm(i)tn(r + jy) ~ Vj cos [20! (r) - 20! (r + jy)] (B12) 
+Vj cos [20i (r) + 20i (r + jy) + Annx - 2-ku] , (B13) 

where we have only retained cosine terms. Written in 
terms of the •§ fields, the first line corresponds to the 
non-Umklapp interaction Eq. (|A15|) with f)W = —Vj. 
The second line is non-oscillatory only at half-filling and 
corresponds to the Umklapp interaction Eq. (|A21[) with 
u {j) = Vj. Finally, terms of the type Eq. (|A"26|) at half- 
filling with q = 1, m = even, arise from Umklapps like 

cos [E£=i20i(r+#)\ 

Let us remark about the effects of compactness of the 
microscopic gauge field. As is known from the ID folk- 
lore, allowing cosines of the dual fields in the action effec- 
tively allows vortices in the microscopic phase variables 
and provides a faithful treatment of the compactness of 
the phase variables. In this respect one may wonder 
about the status of our theory once we allow the de- 
scribed cosine terms in the 0i and 9 2 variables. It turns 
out that it is not complete yet, but becomes so after we 
allow terms like Eq. (|A14[) which in the parton-gauge 
setup correspond to allowing monopoles in the micro- 
scopic gauge field. Since the insertions Eq. (|A14[) have 
ultra-short-range correlations at the EBL fixed pointji 
the issue of monopoles can be safely ignored in the sta- 
ble EBL theory (but of course they cannot be ignored if 
the EBL becomes unstable and the parton fields acquire 
gaps). 

We also note that while we used bosonic partons, a the- 
ory similar to that in Eq. (IB3 | would arise for a so-called 
extremal DLBL state of Ref. |22|, where on the mean field 
level one starts with fermionic partons with flat Fermi 
surfaces (in the context of ring models with K\ < 0). 
Thus, the EBL can be viewed as a solution of such a 
special parton-gauge system, which nevertheless already 
has remarkable properties such as the non-Fermi-liquid 
specific heati C - TTn(l/T). 
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